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ABSTRACT 

We present an analysis of line-of-sight extinction measurements obtained using data from the 
Magellanic Clouds Photometric Survey (Zaritsky, Harris, & Thompson 1997, AJ, 114, 1002), 
which provides 4-filter photometry for millions of stars in the Large Magellanic Cloud. We find 
that visual extinctions are typically larger by several tenths of a magnitude for stars with effective 
temperatures > 12000 K, than for stars with effective temperatures between 5500 K and 6500 
K. Several repercussions of this population-dependent extinction are discussed. In particular, 
LMC distance measurements that utilize old stellar populations, but use extinctions derived 
from OB stars, may be biased low. As a specific example, we show that the LMC distance 
modulus derived from field red clump stars is revised upward relative to published measurements 
by ^ 0.2 mag if one uses the extinction measured for a matched stellar population. Conversely, 
measurements that utilize the youngest stars are subject to greater, and more variable, extinction 
leading preferentially to results that may be biased high. Population-dependent extinction affects 
the interpretation of color-magnitude diagrams and results in an effective absorption law that 
is steeper than that intrinsic to the dust for unresolved stellar systems. We further explore the 
relation between the stellar populations and dust by comparing our extinction map to the lOO/xm 
image of the region and identifying potential heating sources of the dust. We find that although 
regions of high lOO^m flux are associated with young stars, young stars are not necessarily 
associated with regions of high 100/xm flux and that ~ 50% of the 100/xm flux is emitted beyond 
the immediate regions of high OB stellar density. We conclude that 100/im flux should be 
used with caution as a star formation tracer, particularly for studies of star formation within 
galaxies. Finally, we reproduce the observed extinction variation between the hot and cold stellar 
populations with a simple model of the distribution of the stars and dust where the scaleheight of 
the cooler stars is :» than that of the dust (which is twice that of the OB stars; Harris, Zaritsky, 
& Thompson 1997, AJ, 114, 1933). 

Subject headings: dust, extinction — galaxies: distances and redshifts — galaxies; photometry — Mag- 
ellanic Clouds 

1. Introduction HZT), we applied the classic technique of deter- 

mining line-of-sight extinction from the reddening- 
The distribution of dust within galaxies and f^ee colors of OB stars (Johnson & Morgan 1953) 

its effect on our observations are poorly deter- ^o produce an extinction map of a 1.9° x 1.5° re- 

mined and have serious implications for many ar- gjon of the Large Magellanic Cloud (LMC) and to 

eas of extragalactic astronomy. In a previous pa- investigate the relative line-of-sight distributions 

per (Harris, Zaritsky, & Thompson 1997; hereafter of the dust and OB stars. We found that the OB 

stars in the LMC have a scaleheight that is about 
half that of the dust. The combination of this 
result and the observed spatial variations among 
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different stellar populations (cf. Harris and Zarit- 
sky 1999) suggests that each stellar population will 
have its own distinct distribution of extinction val- 
ues. Can the relative distributions of stellar popu- 
lations and dust be mapped and what effects does 
this variation have on observations of the Large 
Magellanic Cloud in particular and of galaxies in 
general? 

The correlation between dust and certain stellar 
populations is well-known {e.g., OB stars and dust 
are preferentially found in galactic spiral arms), 
but generally ignored when treating the global 
effect of dust on galaxy photometry due to the 
complexity and ambiguity such a treatment intro- 
duces. One specific example where such a corre- 
lation may be important is found in the measure- 
ment of the distance to the Magellanic Clouds. 
Different reddenings toward Cepheid, RR Lyrae, 
RGB tip stars, and red clump stars would result 
in apparent differences among the derived distance 
moduli for each of these populations. Apparent 
differences of a few tenths of a magnitude in dis- 
tance moduli are observed (eg. m — M ~ 18.3 for 
RR Lyrae and red clump stars; Layden et al. 1996; 
Udalski et al. 1998; Stanek et al. 1998; Cole 1998, 
Girardi et al. 1998; and m - M ~ 18.5 to 18.7 
for Cepheids; Madore and Freedman (1998) and 
Feast and Catchpole (1997)) and render the dis- 
tance measurements to this key calibration galaxy 
systematically uncertain to about 10%. Although 
there are many factors that give rise to the dis- 
crepancies, the differences, at least in part, may 
arise from applying extinctions derived from a dif- 
ferent stellar populations than that being used as 
the standard candle. 

For unresolved galaxies, even the most detailed 
reddening corrections are fairly rudimentary (typ- 
ically assuming uniform, symmetric distributions 
of dust; cf. Tully & Foque 1985). Although 
increasingly more detailed treatments are being 
constructed and applied (Huizinga 1995; Witt & 
Gordon 1996; Jones, Davies, & Trewhella 1996; 
Kuchinski et al. 1998), these focus primarily on 
the effect of a clumpy absorbing media and have 
not yet included the additional complexity of var- 
ied distributions of dust and stellar populations. 
A complete treatment is obviously a difficult, com- 
plicated problem for which observed examples 
of dust/star geometries would provide invaluable 
constraints. The Clouds provide an opportunity 



to measure the effects of a non-uniform dust dis- 
tribution on different stellar populations through 
line-of-sight extinction measurements for millions 
of stars. 

Once the distributions of stars and dust are 
determined, we can begin to investigate issues 
regarding detailed physics of the interstellar 
medium. A specific example that we can address 
is the outstanding issue of how much of the far 
infrared emission comes from dust heated by the 
UV flux from nearby OB stars (Young et al. 1986; 
Devereux & Young 1990) rather than by the galac- 
tic interstellar radiation field (Lonsdale-Persson & 
Helou 1987, Walterbos & Schwering 1987). The 
resolution of this question impacts the degree to 
which the far-IR flux can be used to measure, or 
at least trace, recent star formation activity. At 
least two dust components are acknowledged to be 
contributing to the far infrared flux: a warm com- 
ponent associated with H II regions (and hence 
massive stars) and a cold component, the cir- 
rus, that is heated by the interstellar radiation 
held. Previous studies have by necessity focused 
on galaxies with marginally resolved stellar pop- 
ulations {e.g., M33; Devereux, Duric, & Scowen 
1997 ) or unresolved populations (cf. Young et 
al. 1996; Lonsdale-Persson & Helou 1987), so it 
has been difflcult to investigate the detailed phys- 
ical connection between stars and dust emission. 
The contribution from the cirrus component is 
thought to depend on Hubble type (Sauvage & 
Thuan 1992) decreasing to a contribution of ~3% 
in Sdm galaxies, like the LMC. The expectation 
that little 100/xm emission will be unassociated 
with young stars in the LMC is simple to test be- 
cause regions populated by young stars are easily 
identiflable from our photometry, the IRAS im- 
age has resolution comparable to the sizes of star 
forming regions, and our extinction map provides 
a measurement of the dust column density. 

With line-of-sight extinction measurements to 
a wide set of stellar populations within the LMC, 
we can begin to explore the relative distribution 
of dust and stellar populations, at least in this 
one galaxy. This study is a necessary step to the 
accurate detailed interpretation of the rich color- 
magnitude diagrams of the Clouds currently be- 
ing produced. To motivate this study in a wider 
context, we have described a few ways in which 
a dust map of the LMC relative to various stel- 
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lar populations may address issues in extragalactic 
astronomy and ISM physics. The previous study 
of extinction in the LMC with the highest spatial 
resolution utilized 2069 O and B main sequence 
stars drawn from a dataset of over 1 million stars 
(HZT). To extend the range of stellar types in- 
cluded in the extinction analysis and enlarge the 
sample, we develop and describe a technique that 
measures the extinction by fitting stellar models 
to the stellar spectral energy distribution. We de- 
rive estimates of both the reddening and effective 
temperature of every star in our survey for which 
4-band photometry is available. In Section 2 we 
briefly describe the data from a region of about 
4° X 2.7° in the LMC. In Section 3 we describe our 
analysis and in Section 4 our results. 

2. Data 

These data come from the ongoing Magellanic 
Clouds Photometric Survey described initially by 
Zaritsky, Harris, & Thompson (1997; hereafter 
ZHT) which is a UBVI imaging survey conducted 
with the Las Campanas Swope (Im) telescope, 
the Great Circle Camera (Zaritsky, Shectman, & 
Bredthauer 1996), and a 2048x2048 pixel CCD 
with pixel scale of 0.7 arcsec pixel" ^. The data 
reduction is described in the original paper except 
for one modification. 

We revise the photometric registration algo- 
rithm because the results were unsatisfactory 
when the original scheme was applied to the 14 
drift scans included here. We developed a superior 
alternate algorithm based on adjusting the pho- 
tometry of the individual subscans reduced with 
DAOPHOT II (Stetson 1987). We assume that 
the photometry of each subscan can be related to 
the true photometry by a single zero point offset, 
and we assume that on average the photometry 
from the survey is accurate (i.e. has a net zero 
point of 0.0 mag). We calculate the relative zero 
point for a particular subscan by evaluating the 
median of the photometric shifts relative to all of 
the neighboring subscans (calculated using stars 
in common within the overlap regions between 
adjacent subscans). Cycling through all of the 
subscans, we identify the subscan with the largest 
median offset. Its photometry is adjusted by that 
offset, and the algorithm begins again. We stop 
iterating when the the largest subscan offset is 



less than 0.02 mag (a value that we judge to be 
below the noise of this procedure). The principal 
remaining difficulty is an artifact of the peculiar 
subset of data currently available. The presence 
of a gap in coverage down the middle of the region 
(cf. Figure 1) renders the relative zero points of 
the right and left half of the region less than ideal 
(and possible issues raised by this deficiency will 
be discussed). The gap will be filled as the survey 
progresses and the entire region described in this 
paper will be surrounded by adjacent scans. 

3. Measuring Extinction 

To measure the extinction we fit three param- 
eters, the effective temperature, the himinosity, 
and the extinction, to the four-filter photometry of 
each star. Because we only have four data values 
to fit per star, more sophisticated modeling (eg. 
an exploration of the extinction curve) is unwar- 
ranted. The requirement that a star be detected 
in all four bands eliminates a substantial fraction 
(0.60) of the entire sample of stars, because the U 
and I band images are shallower than the B and V 
images. Even though the procedure appears rela- 
tively straightforward there are various subtleties 
that must be addressed. 

First, we must adopt a set of suitable stellar 
models to fit. We use the library of models pro- 
vided by Lejeunne, Cuisinier, and Buser (1997). 
These models span a wide range of metallicitics, 
effective temperature (Tg, and loosely referred to 
as the temperature throughout this paper), and 
log g. Special care was taken by those authors 
to place the models on the correct color scale so 
that the theoretical solar-abundance model yields 
the observed solar U through L colors. This cor- 
rection is particularly important for a use such as 
that being described here. Of their many libraries, 
we use the solar metallicity library (IcbpOO.cor) for 
calibration and for a test of the metallicity depen- 
dence of the method, the [m/H] = —0.5 library 
(lcbm05.cor) for the LMC analysis (the average 
[0/H] gas phase abundance among LMC HII re- 
gions is —0.53; Pagel et al. 1978), and the [m/H] 
= — 1 library (ImbmlO.cor) for a test of the metal- 
licity dependence of the algorithm. 

Second, we must adopt a conversion system be- 
tween magnitudes and fluxes. One option is to 
adopt the absolute flux densities of an AOV star 
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from an existing reference, such as Straizys (1992). 
In principle, the standard star's flux densities en- 
able us to convert any flux into a magnitude, or 
vice versa. However, our photometric system may 
not be an exact match to this particular photomet- 
ric system or to that adopted for the theoretical 
models. Instead, we adopt a more self-consistent 
approach by defining one of the models to repre- 
sent an AOV star {Te = 9500° and log(s') = 4). 
We calculate fluxes through the our filter system 
and assign these fluxes to a mag star in each fil- 
ter. With this calibration, fluxes and magnitudes 
through each of the four passbands can be evalu- 
ated for any temperature-surface gravity model. 

Third, we must adopt a function for the wave- 
length dependence of extinction across the optical 
portion of the spectrum. Our adopted extinction 
law is a standard Galactic curve (Schild 1977 and 
references therein), which for A > 2600A is vir- 
tually identical to the LMC extinction curve even 
in environments that show signiflcant differences 
blueward of 2600A (Nandy et al. 1980; Nandy, 
Morgan, & Houziaux 1984). We parameterize the 
normalization of the extinction law using the ex- 
tinction in the visual band. Ay. Variations in Ry, 
Av/E{B — V), are known to occur in the galaxy 
(cf. Carrasco, Strom, & Strom 1973; Cardelli, 
Clayton, & Mathis 1989), but we cannot resolve 
such variations in the LMC with our data. Koorn- 
neef (1982) demonstrated that the typical value of 
Rv for the LMC is similar to that of the typi- 
cal line of sight in the Galaxy, 3.1 (which is the 
adopted value for this study). 

The temperature and line-of-sight extinction of 

a particular star arc then derived by identifying 
the best fit model for that star. The algorithm 
compares the observed flux in each passband to 
that calculated from a model. The model is renor- 
malized to minimize residuals and is calculated 
using the observational uncertainties. We allow 
the flux to be rcnormalizcd to avoid assuming a 
distance to the star (which is a more critical free 
parameter for the local calibrating sample to be 
discussed than for the LMC stars). The residu- 
als are calculated as fractional differences between 
the observed and predicted flux {i.e. Amag, which 
weighs each filter in the determination of more 
equally than if Afiux is used). We cycle through 
all available stellar models with T > 3000 K and 
for each model, we cycle through < Ay < 2 in 



steps of 0.01 mag. The parameters of the model 
with the lowest value of a^re accepted as the 
stellar temperature and Ay. 

We use local stars of known spectral type 
and magnitude to calibrate and test the proce- 
dure. Stellar data, UBVI magnitudes and spec- 
tral types, are collected from a series of catalogs 
available on-Une at the Strasbourg data center 
(http://cdsarc.u-strasbg.fr/; Lanz 1986; Mermil- 
liod 1986; Houk & Smith-Moore 1988) for 107 
HD stars. Using the relationship between stellar 
type and effective temperature given by Straizys 
(1982), we convert the spectral type into an ef- 
fective temperature for each star. In Figure 2 we 
show the results from various tests of fitting the 
stellar models to these data. In the upper left 
panel of the Figure, we have plotted the derived 
effective temperatures versus the known tempera- 
tures, which we refer to as spectroscopic tempera- 
tures, obtained by assuming zero extinction along 
the line of sight and no photometric offset between 
the theoretical photometric system of the models 
plus our filter transmission curves and that used 
to produce the published magnitudes. 

The agreement between the two temperatures is 
encouraging but there arc two qualitative failures 
that must be corrected. First, the derived tem- 
peratures are systematically lower than the spec- 
troscopic temperatures. Second, the deviations 
become large for stars with Te > 10000 K. In 
the upper right panel, we plot the results allowing 
the algorithm to fit extinctions. The agreement is 
significantly better. There is no noticeable offset 
between photometric and spectroscopic tempera- 
tures for stars with Te < 10000 K and the situa- 
tion has improved somewhat for the hotter stars 
as well. Next we allow the photometric zero points 
of each filter to change, and we select the best off- 
sets by minimizing the residuals about the 1:1 line. 
The results obtained with these offsets are shown 
in the lower left panel. We have improved the fit 
at the hotter effective temperatures, although this 
improvement has been obtained by construction. 
The photometric offsets are modest (1% in U and 
B, 7% in V and 12% in /). These offsets may in 
part arise from a mismatch of the reference model 
and an AOV star or from differences in the filter 
transmission curves between photometric systems. 
Finally in the last panel (lower right), we plot the 
result of applying the algorithm (with freedom to 
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fit extinction and with the new photometric off- 
sets) to data for UBVRI standards taken from 
the Astronomical Almanac (1981). Again, we have 
converted from stellar types to Te using the con- 
version table from Straizys (1982). The photomet- 
ric temperature agrees well with the spectroscopic 
temperature over the range 3000 K to 35000 K. 
The agreement between the fitted and published 
temperatures for this independent sample demon- 
strates that the photometric offsets derived using 
the first sample are generally appropriate for other 
data as well. We conclude that we can recover Te 
to a fractional precision of ^ 10% over this tem- 
perature range. 

Although we obtain estimates of the line-of- 
sight extinction toward these standard stars, there 
are no published measurements to compare with. 
Instead, we test our extinction measurements by 
artificially reddening the photometry of these stars 
and attempting to recover the added extinction. 
This approach tests the degeneracy among vari- 
ous stellar models plus reddening, but does not 
test the extinction law itself. We apply the same 
supplemental visual extinction (0.2, 0.4, and 0.8 
mag) to the 107 standard stars, refit the mod- 
els, and accept the best fits only if < 3 (as 
we will do with the analysis of LMC stars). In 
Figure 3 we plot a running median of the newly 
derived extinction minus the initial derived extinc- 
tion (bins of 11 stars, with stars sorted by photo- 
metric temperature;) and compare the result with 
the input supplemental extinction, which is plot- 
ted as a dotted line. For low supplemental ex- 
tinction values. Ay = 0.2, which corresponds to 
low total extinction because the initial extinctions 
are generally small, the supplemental extinction is 
recovered well across the entire range of tempera- 
tures. As the extinction increases, we notice some 
serious regions of failure. In particular, the region 
of photometric Te between 7000 and 9000 K has 
extinction values that are significantly underesti- 
mated for supplemental extinctions of 0.4 and 0.8 
mag. 

To examine the failure of the method for stars 
of intermediate temperatures, we plot the results 
for supplemental extinction of 0.8, but sort stars 
by spectroscopic, rather than photometric, tem- 
perature (lowest panel of Figure 3). By compar- 
ing the bottom two panels, in particular the tem- 
perature range over which the Ay determination 



fails, we can understand what has happened. Stars 
with large extinctions and Te ~ 10000 K are well 
fit by models with cooler stellar atmospheres and 
smaller Ay- This problem is a manifestation of 
the classic degeneracy evident when applying the 
reddening- free method in the U — B,B — V color- 
color diagram. For small extinctions, the degen- 
eracy is avoided, but for large extinctions various 
combinations of Te and Ay are possible for stars 
with true Tg's between 7000 and 12000 K. We also 
see a failure mode below 5000 K for our algorithm. 
We suspect that the problem at these low tempera- 
tures is a combination of insensitivity to reddening 
in stars that are already quite red and the effect 
of strong molecular absorption features that effec- 
tively introduce noise into the process of matching 
the models to the observations. We conclude that 
for the range of extinctions expected in the LMC, 
we should constrain our analysis to stars with pho- 
tometrically derived Tg's between 5500 and 6500 
K and greater than 12000 K. 

We run the algorithm on the stellar catalog 
obtained for our survey region using the [m/H] 
= —0.5 models and the adopted photometric off- 
sets. We only use stars with moderate or better 
photometric errors (cry < 0.2, as < 0.1, ay < 0.1, 
and a I < 0.1) and apply magnitude cuts (exclude 
stars with U and B > 22 and F and / > 21). We 
obtain satisfactory fits (x^ < 3) for 94% of the 
stars that satisfy our criteria (a total of 543,860 
stars within this region). For our subsequent anal- 
ysis we exclude stars with fitted Te < 5500 K and 
6500 K < Te < 12000 K. In the lower allowed 
temperature range we have 39,613 stars and in the 
upper allowed temperature range we have 106,181 
stars. 

The various photometric selection criteria pose 
a potential risk of introducing unwanted biases. A 
limiting magnitude necessarily implies that heav- 
ily extincted stars will be lost from the sample. 
The selection based on photometric errors simi- 
larly introduces a "fuzzy" magnitude limit because 
the errors correlate strongly with magnitude. Ex- 
cept for extremely extincted stars, these limits are 
unlikely to be important for the hotter stars be- 
cause those stars are at least a magnitude brighter 
than the V band limit and have B — V ^ 0. There 
is greater potential for biases in the colder sam- 
ple because these stars are fainter and redder, and 
so more likely to be lost from the U band data. 
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which is our shallowest. Upon inspection we find 
that although the magnitude limit {U > 22) is well 
below (~ 2 mag) the red clump population, which 
is where the bulk of the cold stars lie, the photo- 
metric error cut effectively raises the limiting mag- 
nitude to ~ 21 mag. To test whether this seriously 
affects the recovered distribution of extinctions we 
have selected a second sample of cold stars with 
0.2 < au < 0.4 to recover some of the excluded 
stars. The comparison of extinction values among 
the two samples is presented in Figure 4. There 
is no indication that high extinction (> 0.4) stars 
are being artificially dropped from the sample due 
to the au < 0.2 cut. We conclude that although 
there is no way to avoid losing the fainter stars 
along lines-of-sight with extreme extinctions (> 1 
mag), selection criteria are not significantly bias- 
ing the recovered distribution of Ay^s. 

Using the derived line-of-sight extinction val- 
ues, we construct an extinction map by spatially 
grouping the measurements. We regrid the sur- 
vey region and begin with a map pixel that cor- 
responds to 120 original pixels (84 arcsec). Re- 
gions of high stellar density could be treated with 
smaller map pixels, but this choice of pixel size 
was found to be a good compromise for the en- 
tire survey region. For map pixels that contain 
at least three stars with extinction measurements 
the median extinction value is adopted. For map 
pixels with two or fewer values, we increase the 
size of the map pixel size by another 120 original 
pixels. If fewer than three stars are found within 
this larger map pixel, we increase the map pixel 
size and try again. We continue incrementing the 
map pixel size until the new, larger pixel contains 
at least three extinction values. All original pix- 
els contained within the enlarged map pixels are 
assigned the median value. This algorithm results 
in a final map that has regions of high and low 
spatial resolution, but relatively uniform signal- 
to-noise. Finally, we create a mask that sets any 
original map pixel containing zero usable stars to 
zero extinction because such pixels are likely to be 
either off the surveyed region or to be corrupted 
by technical problems. 

The two maps resulting from the analysis of 
high and low Te stars are presented in Figure 5 
(with higher extinction regions shaded darker). It 
is evident that the average extinction values are 
larger for the hotter stars. We will return to 



this point in Section 4. We have marked, with 
a vertical line to the right of the images, the scan 
that is known to have 7-band data taken in non- 
photometric conditions and we suspect that the 
feature at the far right of that scan (near the right 
edge) is a result of these non-photometric data. 
These non-photometric data, through the process 
of photometric registration discussed in §2, may 
also have partially corrupted the data adjacent to 
the top of the non-photometric scan, which ap- 
pear to have an anomalously low Ay in the up- 
per panel of the Figure. Finally, the other feature 
we note in our preliminary discussion of the maps 
is the left/right offset present in the lower panel. 
We cautioned previously that the central coverage 
gap renders the relative photometric normaliza- 
tion of the left and right sections somewhat un- 
certain (0.05 mag) and a slight shift in the mean 
extinction is seen between the two sides for the 
cool stars. Other than these issues, there is some 
slight structure that correlates with subscan or 
scan edges. 

4. Discussion 

The extinction maps shown in Figure 5, in 
particular the one derived from the stars with 
Te > 12000 K, show large coherent, filamentary 
structures. Most noticeable is the broken ring of 
extinction in the left hand of the survey region that 
surrounds the central stellar association of LMC-4. 
Most of the high extinction regions visible in the 
top panel of Figure 5 can also be seen as regions 
of obscuration in the stellar density map (Figure 
1). This coincidence confirms the features in the 
extinction map as real regions of enhanced extinc- 
tion. Because of the larger number of stars used, 
we can trace features at finer spatial resolution 
than done by HZT. After examining the sensitiv- 
ity of the algorithm to the choice of mctallicity, 
we proceed by comparing the reddening map to 
previous reddening maps and other observations 
of the same region. 

4.1. Variation with Metallicity 

Although the difference in the inferred extinc- 
tion values between high and low Te stars is sig- 
nificant, we must exclude the possibility that this 
difference is due to physical reasons other than 
a difference in line-of-sight extinction. One such 
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possible reason is the metallicity difFcrcncc that 
presumably exists between the younger and older 
stellar populations. We have selected to fit mod- 
els that correspond to the current gas-phase abun- 
dance in the LMC for all of the stars, which may 
not be appropriate for the low Te stars. We test 
the technique's sensitivity to mismatched mctallic- 
ities between the stars and models by reapplying 
the algorithm using tenth-solar and solar metallic- 
ity models on a random subset of the data. 

The resulting differences in and Ay as a 
function of Te are shown in Figure 6 (we plot the 
difference in the inferred quantity between a run 
using the standard [m/H] = — 0.5 models and the 
other [m/H] models). There are no systematic 
differences in either the effective temperature or 
extinction, except possibly at temperatures below 
5000 K for Av when using the solar metallicity 
models. In both plots of the derived extinction 
differences, there are one or more points, within 
our Te range, for which the average difference 
is large (> 0.1 mag). For those points, we have 
also evaluated the median value within those bins 
and plotted these values as open circles. The me- 
dian, which we use when generating our extinction 
maps, is consistent with zero difference. A change 
of a factor of 10 in the adopted abundance appears 
to have little or no net effect on the derivation of 
Av and Te over large samples of stars. This result 
reflects some of the difficulty in obtaining metallic- 
ity measurements for individual stars from broad 
band colors, and assures us that metallicity mis- 
matches, which are certainly present at some level, 
are not responsible for the offset seen in Ay be- 
tween high and low Te stars. 

4.2. Comparison to HZT's OB Star Ex- 
tinction Measurements 

A straightforward test of the derived extinction 
values for the high Te stars is a direct compar- 
ison with previous results. HZT compared their 
results to previous studies (primarily. Hill et al. 
1994, Masscy et al. 1995, Oestreicher & Schmidt- 
Kaler 1996) and found that if similar selection cri- 
teria are used the distribution of reddening values 
agree. Because the HZT results are based on a 
subset of the data used here, the comparison to 
the HZT results that we describe below directly 
tests our algorithm. The agreement between the 
HZT results and previous studies address the reli- 



ability of the data. 

We match the 1615 stars in common between 
the two studies and present the comparison in Fig- 
ure 7. On average the two measurements agree 
well and there is no indication of a bias between 
the two techniques, but the dispersion about the 
1:1 line is large, 0.21 mag. Some of this scatter 
is attributable to the different photometric zero 
pointing technique used here and correlates with 
position in the survey region. The scatter among 
extinction mcasiirements within smaller areas of 
the survey is about half of the global scatter and 
this reduced local scatter gives rise to the appear- 
ance of two tighter sequences of stars above and 
below the 1:1 line. From this comparison, we con- 
clude that the uncertainty in Ay for a single star is 
at least as large as 0.1 mag and possibly as large as 
0.2 mag. The comparison demonstrates that the 
reddening to any individual star is significantly un- 
certain and the results from these techniques must 
be averaged over several stars, as we do to produce 
our extinction maps, to mitigate these uncertain- 
ties. Nevertheless, we believe that the newer pho- 
tometric zero-pointing technique is superior and 
so the current measures of Ay are preferable over 
the HZT values. For the averaging process to be 
physically meaningful, the region over which the 
averaging is done must have homogeneous extinc- 
tion properties and the final extinction map will 
not resolve structures on scales smaller than the 
size of the region chosen for averaging or median- 
ing. 

4.3. Comparison to Published Extinctions 
for OB Regions 

We compare our reddening values with those 
published for OB regions within our survey. We 
have identified four Lucke-Hodge regions ((LH 38, 
54, 83, and 114; Lucke and Hodge 1970) in com- 
mon between our survey and the studies by Hill, 
Madore, and Freedman (1994) and Oey (1996). 
Our median values of E{B — V) for stars in the 
four regions are 0.12, 0.13, 0.13, 0.14, respectively. 
The published values of E{B - V) are 0.15 (Oey), 
0.01 (HiU et al.), 0.16 or 0.11 (Hill et al. and Oey, 
respectively for LH 83), and 0.08 (Oey) for the 
four regions in the same order as above. The pub- 
lished values are typically based on only tens of 
observed stars and have large uncertainties (Hill et 
al. quote uncertainties of ^ 0.08). The values we 
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attribute to Ocy arc medians taken from her list of 
reddenings to individual stars in the regions. We 
find no systematic differences between our values 
and the published ones, although the difference in 
one case (LH 54) is slightly larger than the uncer- 
tainty (but only by a factor of 1.5). This differ- 
ence may be unimportant because Hill et al. note 
that their data for this region contained signifi- 
cant scatter among reddenings and many negative 
reddening values, suggesting that there may have 
been a problem with the data. Other than the LS 
54 field, the others are entirely in agreement with 
the published data, although the large uncertain- 
ties preclude more detailed comparisons. 

4.4. Comparison to the OGS's Foreground 
Extinction Measurements 

We proceed to test whether the extinctions de- 
rived for low Te stars are consistent with previous 
studies. The only comparable study is that by 
Oestreicher, Gochermann, & Schmidt-Kaler 1995 
(hereafter, OGS) which used Galactic stars in the 
direction of the LMC to measure the foreground 
extinction toward the LMC. We divide our low Te 
sample into foreground Galactic stars and LMC 
stars by attributing the large vertical plume of 
stars seen in the Hess diagram at i? — ^ ~ 0.6 and 
y < 18 to Galactic stars (cf. Figure 1 from HZT, 
and references therein) . By examining the extinc- 
tion values of the foreground and LMC popula- 
tions separately and comparing them to the mea- 
sured foregroimd (Galactic) extinction from OGS, 
we can determine whether (1) our foreground val- 
ues are in agreement with OGS, (2) whether our 
LMC values are consistent with (i.e., >) the fore- 
ground extinction values, and (3) whether the 
metallicity difference between the Galactic and 
LMC populations grossly affects the derived ex- 
tinctions. 

We begin by comparing our extinction values 
for the Galactic and LMC low Te populations with 
the OGS values. Unfortunately, the comparison is 
not ideal because the OGS data do not sample the 
region as densely as do our data. Figure 8 shows 
the comparison between reddenings to individual 
OGS foreground stars and values inferred from an 
extinction map constructed using our foreground 
stars. There is no correlation visible (although the 
range of values along the axes are similar). The 
scatter demonstrates the large uncertainties in the 



values along any single line-of-sight. 

Although the scatter in individual lines-of-sight 
is large, over large areas the extinction map may 
be reproducible and reliable. To test this possi- 
bility, we first compare the distribution of OGS 
extinction values, including all OGS stars across 
the entire LMC to improve the statistics (no de- 
tectable difference was seen when only using the 
foreground stars within our survey region), to the 
distribution of extinction values derived for the 
foreground stars in our sample in Figure 9. Al- 
though the Ay values from the OGS data are 
somewhat more widely distributed than the val- 
ues we infer from our foreground population, the 
modes of the two populations are similar (both ~ 
0.11 mag). The wider OGS distribution may im- 
ply either that the uncertainties are larger in the 
OGS data than in our data, that we have underes- 
timated Ay for some of the stars in our survey, or 
that the region sampled by OGS included regions 
of higher foreground extinctions. The quoted un- 
certainties in the OGS data are mostly < 0.06 
mag, insufficient to account for the larger spread 
in values. This conclusion is further supported by 
the rough similarity in the fraction of points with 
Av < in both datasets. In bulk, the statistics 
of the foreground extinction as derived from our 
data and OGS's are quite similar. 

We compare the Ay distributions for the Galac- 
tic and LMC stars in the inset of Figure 9. There is 
little difference in the distributions except for the 
more significant tail of higher extinction values in 
the LMC stars. This agreement demonstrates that 
there are no gross differences between the extinc- 
tions derived from the different set of stars due 
to physical differences, like metallicity, and that 
the LMC data are consistent with the extinction 
inferred from the foreground data (from both the 
OGS data and our own). 

Although statistically the extinctions inferred 
from the OGS data and our own are similar, we 
need to confirm whether we can accurately trace 
the spatial distribution of dust to validate out ex- 
tinction map. To improve the statistics for the 
comparison between the OGS observations and 
ours, we compare the extinction map produced 
from the OGS data to that obtained using all of 
our low Te stars (the LMC stars should trace the 
foreground extinction pattern in addition to any 
pattern imposed by dust internal to the LMC). 



8 



Pixel values in our low Te extinction map are 
tagged according to the extinction value given by 
the OGS map for the corresponding location. We 
then bin the data using the OGS extinction val- 
ues and set the bin size to include extinction values 
from at least 200 different pixels within each bin. 
The average of our extinction values within each 
bin is plotted versus the OGS extinction value in 
Figure 10. There is a clear correlation between 
our low Te extinction measurement and the OGS 
data, demonstrating that the bulk properties of 
the OGS map and our extinction map are repro- 
duced and that our analysis is recovering the spa- 
tial behavior of foreground extinction in this re- 
gion. There are a few discrepant values in bins at 
low values of the OGS extinction, but the pixels re- 
sponsible are localized in our map. Because of this 
spatial localization, the discrepancy is either due 
to slightly erroneous photometry in a limited re- 
gion {e.g., poor aperture correction) or to a patch 
of extinction internal to the LMC extinction that 
is not represented in the OGS foreground extinc- 
tion map. We conclude that the strong correlation 
demonstrates that we recover the spatial distribu- 
tion of extinction and that there is no dominant 
component of extinction internal to the LMC that 
affects the low Te stars. 

Despite the strong correlation, the relation be- 
tween the OGS results and ours is similar to that 
seen in Figure 9 (our extinction values < OGS 
extinction values). This offset may be the result 
of a slight bias in our method when applied to 
the cool stars {e.g., inclusion of some intrinsically 
higher Te stars for which Ay is underestimated or 
the exclusion of some higher Ay stars due to our 
photometric selection), or it may be indicative of 
unresolved structure in the OGS map. If areas of 
high extinction are small on the sky, then a map 
produced from sparsely sampled data would oc- 
casionally include large inferred areas of high ex- 
tinction (large areas of low extinction would also 
be present but would be a better representation 
of "reality" ) . When such a map is compared to a 
map with higher resolution, a bias in the observed 
sense would exist. We cannot determine which 
of these two scenarios is responsible for the sys- 
tematic discrepancy we observe, but we conclude 
that to within this possible slight bias our extinc- 
tion determinations for low Te stars recover the 
statistical properties and spatial behavior of the 



foreground extinction found in previous determi- 
nations of the foreground extinction. 

4.5. Comparison to the IRAS lOO/im Im- 
age 

To determine whether the medianing we per- 
form to mitigate the uncertainties and create the 
extinction map results in a physically meaningful 
map, we compare our map to images of the far in- 
frared emission. The far IR emission comes from 
the dust, but the emitted flux depends on the dust 
temperature and hence on the local radiation field. 
Therefore, although lOO/um flux is expected to be 
a tracer of dust, we do not expect it to be a com- 
plete map of the extinction. We have extracted the 
IRAS lOOjum image of the area and present it, in 
conjunction with the extinction map from the high 
Te stars and the density-weighted image of main 
sequence stars with V < 16.5 and B — V < 0.5, 
in Figure 11 (B3 and earlier main sequence stars). 
We apply our mask to the IRAS image to produce 
an image that exactly matches the area for which 
we have optical data. Most of the prominent dust 
features are visible in both the extinction map and 
the IRAS data, which further confirms that struc- 
tures seen in the extinction map match true dust 
structures. 

There are only a few instances where signifi- 
cant 100//m flux exists without strong extinction 
or where strong extinction exists without signifi- 
cant lOO/Lim flux. Most of the differences can be 
traced to the presence or lack of nearby young 
stars. For example, the extinction map shows a 
long vertical filament in the upper portion of the 
left half of the survey region. In extinction, this 
filament does not have strong variations along the 
filament. However, the lOO/im emission is strongly 
peaked near the lower portion of the filament, co- 
incident with the presence of young stars that pro- 
vide the UV flux needed to heat the dust. A 
more extreme example is the small, brightest knot 
of lOO/Ltm emission in the upper half of the left 
side of the region, which has no noticeable extinc- 
tion consequence but is coincident with a group of 
young stars. Highly localized dust may be diluted 
in our map by the medianing we do, but because 
this area includes a large number of OB stars we 
have not diluted the extinction at this specific lo- 
cation. These examples demonstrate that both 
the column density of dust and the presence of 
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heating sources are important in determining the 
lOO/um morphology, but that in general the ex- 
tinction features determined from medianing over 
several stars per pixel are confirmed by the 100/um 
emission. 

Although strong IR emission appears to re- 
quire the presence of coincident young stars, young 
stars are not always indicative of a strong lOO/im 
emission feature. The dominant concentration of 
young stars (of both the O and B stars) within the 
survey region is in the middle of the left side of the 
region. However, this area is a one of the lowest 
areas of 100/zm flux, presumably because there is 
less dust there than in other areas — as suggested 
by the Ay map. Perhaps the young stars have 
removed the dust along the line of sight by blow- 
ing a bubble out of the disk plane in our direction 
(cf. MacLow and McCray 1988). Such a scenario 
would also be in agreement with the ring-like dust 
morphology surrounding this region seen both in 
the extinction map and the 100/zm image. Obser- 
vations of the ionized gas in this region are also 
qualitatively consistent with this scenario (Hunter 
1994). If the bubble is as large in the line-of-sight 
direction as in the cross direction (~ 1 kpc), then 
it has blown out of the LMC disk plane. 

We quantify the connection between dust, 
lOO/um flux, and young stars by examining the 
correlation among pixel values in the various im- 
ages. By identifying pixels in the lower panel of 
Figure 11 that contain 2 or more OB stars (about 
16% of all pixels), we find that 71% of all OB 
stars are contained in these pixels, but that only 
21% of the lOO/um flux is contained within these 
same pixels. Each pixel corresponds to a projected 
size of 20 pc (the dust map is smoothed using a 
Gaussian of 3 pixels full- width at half maximum). 
Although there is a slight increase in 100/im flux 
over that expected from a uniform distribution 
(21% vs. 16%), most of the IR flux is not coming 
from the regions of highest OB star concentra- 
tions. Pixels without a single OB star contain 
54% of the lOO/xm flux. If instead we divide the 
region on the basis of pixels that have higher than 
median 100/;im emission (12.1 MJy/sr; 50% of the 
pixels), we flnd that this area contains 75% of 
the OB stars. These numbers quantify the visual 
impression from Figure 11 that regions with the 
highest lOO/xm flux lie near sites of OB stars, but 
that the local presence of OB stars does not im- 



ply strong 100/im emission. These arguments are 
complicated by the unknown three-dimensional 
nature of the problemk and the fact that OB stars 
can heat dust hundreds of pc away. 

The relative distribution of young stars and 
100/xm flux within this region illustrates why it 
has been difficult to resolve the issue of the heat- 
ing source of dust in other galaxies. Although the 
lOO/Ltm flux correlates in general with young stars 
[i.e., more 100/im flux implies more young stars), 
a substantial fraction of the 100/im flux is diffuse 
(i.e., coming from beyond the sites of OB stars) 
and there are regions rich in OB stars with little 
or no 100/im flux. It is evident from the distribu- 
tion of the two components that 100/im flux can- 
not be used to trace star formation in detail (even 
if globally correlated), and that there may be a 
large dispersion in global characteristics as well, 
although we cannot address that issue with these 
data. Our data directly demonstrate the scatter 
in local properties, but global properties may still 
be correlated if the galactic radiation field is dom- 
inated by the radiation from OB stars. The local 
scatter in any relation between young stars and 
100/im flux is evident if one considers the central 
region of LMC 4 (which is currently nearly empty 
of dust and hence of 100/im emission). The main 
arc of young stars in this region has 17 pixels that 
contain as many OB stars as the densest pixel cor- 
responding to the bright 100/im emission knot just 
above the arc. If there was dust within the arc of 
young stars one would expect the flux to be at 
least 17 times greater than that coming from the 
knot, which is evidently not the case. Therefore, 
although large 100/im fluxes are likely to be an ac- 
curate indication of star formation, the converse is 
less likely to be true and 100/im flux has a large 
dispersion as a local tracer of recent star forma- 
tion within galaxies. A complete treatment (UV 
photometry to measure ionizing flux, optical, HI, 
and a variety of infrared bands) is necessary to 
make further progress in resolving the full energy 
budget of this process, but is beyond the scope of 
the current study. 

4.6. Variation with Stellar Population 

We have demonstrated that the extinction mea- 
surements for the low and high Te stars are in 
agreement with previous measurements of extinc- 
tion and emission from dust. Therefore, we can 
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begin to compare with some confidence the rela- 
tive distributions of dust that affect the low and 
high Te stars. The upper and lower panels of Fig- 
ure 5 demonstrate that the line-of-sight extinc- 
tions are generally much lower for the low Te 
stars. This impression is confirmed by the Ay 
histograms for the two stellar populations (Fig- 
ure 12). The trend of declining mean extinction 
with declining Te (and presumably with increas- 
ing mean age) is also present within subsets of the 
high Te data. In Figure 13 we plot the Ay distri- 
bution for four subsets of high Te stars. There 
are several trends evident as Te increases: (1) 
the number of very low Ay (< 0.1 mag) stars de- 
creases, (2) the peak of the distribution narrows 
and stays at an almost constant Ay, and (3) the 
tail of high Ay values becomes increasingly more 
significant as Te increases. Even within the set of 
high Te stars, higher values of extinction are more 
common for the hotter, younger stars. To test this 
result, we also plot the HZT results using the clas- 
sical reddening- free method (the Q method). Be- 
cause there are fewer stars in the HZT sample, we 
divide it into only two subsamples at F = 15. As 
hmiinosity (and Te) increase, this sample repro- 
duces both the smaller relative number of low Ay 
values and the larger relative number of high Ay 
values. 

Even though the non-Gaussian Ay tail affects a 
small fraction of the total stellar population, it can 
have serious effects on certain obscirvations and 
suggests that the dust is highly localized. Clumpy 
regions of high extinction effectively flatten the 
extinction curve (in the extreme case where one 
has only completely optically thick clumps, some 
stars will be extincted out of the sample, but the 
colors of the observed stars may not change; cf. 
Szomoru & Guhathakurta 1999). This flattening 
of the effective extinction curve for the LMC as 
a whole is countered by our observations that the 
bluer stars are more heavily extincted than the 
red stars, which steepens the effective extinction 
curve. The degree to which these two effects cancel 
depends on the detailed distribution of dust and 
stars, which are still poorly determined in detail. 
Lastly, because of the non-Gaussian Ay tail, gen- 
eral properties of the Ay distribution, such as the 
mode, may not truly reflect the amount of light 
lost from the entire population due to the dust. 

To demonstrate the error introduced by using 



the mode of the Ay distribution to correct pho- 
tometry, we compare the amount of light extincted 
in the 12000 K < < 14000 K sample to that 
extincted in the 22000 K < < 45000 K sam- 
ple. One might at first believe that fractionally 
the extincted light would be comparable because 
the modes of the two distributions arc quite sim- 
ilar (cf. Figure 13), particularly if the available 
data were insufficient to significantly populate the 
high Ay tail in the hotter sample. However, us- 
ing the full distribution of Ay^s and assuming that 
all of the stars within the temperature range have 
the same luminosity [i.e., no Ay — Te correlation 
within the bins), we find that 29% of the V-band 
light is extincted from the cooler sample vs. 42% 
from the hotter sample. As stated before, this pro- 
gression makes the effective extinction curve for 
the galaxy as a whole (ignoring scattering in un- 
resolved galaxies) steeper {i.e., smaller Ry). We 
conclude that even among high Te stars it is criti- 
cal to measure the distribution of Ay to accurately 
correct for extinction. 

Finally, we compare the extinctions from the 
hot and cool stellar populations by comparing 
the extinction maps pixel by pixel (Figure 14). 
Again, we bin and average (so that each bin con- 
tains values from at least 200 pixels). The posi- 
tive correlation in Figure 14 demonstrates that the 
dust affecting the hot population also affects the 
cool population, but the correlation is not one-to- 
one. To understand this correlation, we examine 
a model of the relative distributions of stars and 
dust in which the LMC dust is in a plane of scale- 
height ^ scalehcight of the cool stars. In such a 
geometry, half of the cool stars are only extincted 
by Galactic foreground dust. The LMC dust par- 
tially fills the central LMC disk plane, so a yet- 
unspecified fraction of the stars behind the LMC 
midplane are extincted by the Galactic foreground 
plus the LMC disk and the remainder arc seen 
through holes in the LMC disk and therefore only 
extincted by the Galactic foreground dust. The 
extinction due to the LMC disk is taken to be the 
extinction predicted from the high Te stars minus 
the mean foreground extinction and constrained to 
be > 0. The foreground extinction is taken to be 
0.173 mag (the average from the extinction map 
constructed from the OGS data for the relevant 
region of the LMC), and we plot in Figure 14 the 
results of adopting dust "filling factors" for the 
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LMC dust of 0.0, 0.5, or 1.0. A filling factor of 1 
corresponds to the case where the disk extinction 
is exactly that given by the high Te extinction 
map, while a filling factor of 0.5 corresponds to 
the case where half of the lines-of-sight are clear 
and half are extincted as described by the high Te 
extinction map. This simple model with a filling 
factor of 1.0 successfully reproduces the behavior 
in the Figure. The extinction predicted by the 
model toward the left of the Figure is lower than 
observed, but only by 0.02 mag. If one adds 0.02 
mag to the foreground model (to resolve this dis- 
agreement), then the models with a filling factor 
~ 0.5 work equally well. In either case, we con- 
clude that a model where the low Te stars have a 
significantly larger scalchcight than both the OB 
stars and the dust and where the dust structure 
recovered from the high Te stars is taken to rep- 
resent that of the midplanc sheet works well in 
reproducing both the low and high Te extinction 
behavior. 

4-6.1. Some Implications for the Distance Scale 

One area where accurate extinction corrections 
are critical is distance scale work. The population- 
dependent extinction observed in the LMC is also 
presumably present in other galaxies and can influ- 
ence observations of those galaxies. As such, it be- 
comes important not only to measure the extinc- 
tion local to the source of interest {e.g., Cepheid or 
supernovae), but also to measure the distribution 
of extinctions for a correspondingly similar stellar 
population. Both of these goals obviously become 
increasingly difficult the more distant a galaxy is, 
but they can even be important for studies of 
the LMC. Inappropriate reddenings adopted us- 
ing mismatched stellar populations can affect the 
derived distance to the LMC, which is often a key 
calibrating galaxy. Reddening estimates are gen- 
erally determined from OB stars (cf. HZT) which 
are bright and for which unambiguous reddening 
measurements can be obtained using the straight- 
forward application of the classical reddening free 
color technique (the Q method). However, the 
results in this paper illuminate several key diffi- 
culties with this approach. First, the reddenings 
determined in this way should not be applied to 
Pop II standard candles (cf. red clump stars and 
red giant tip stars). Second, even among the high 
Te stars, there can be differences in the distribu- 



tion of extinction values. Third, because the dis- 
tribution can be highly non-Gaussian, significant 
errors can be introduced by using certain statisti- 
cal properties (e.g., mode) rather than the entire 
Ay distribution to correct the photometry. 

We do not intend here to provide a complete 
review of how variable extinction will affect each 
distance method. However, a few examples will 
demonstrate that the magnitude of the effect is 
comparable to the degree of disagreement between 
various published results. First we consider the 
use of the LMC to zero-point the Cepheid P-L re- 
lation. We have already demonstrated in Figure 
13 that even among stars with Te > 12000 K (the 
young stellar population and therefore appropri- 
ate for Cepheids) there are differences in the dis- 
tribution of extinction values. Any determination 
of the extinction will be dominated by the coolest 
stars included because there will be many more 
of them. For example, the mean visual extinction 
obtained using all of the stars with Te > 12000 
K is 0.42 and the mean using only those stars 
with Te > 22200 K is 0.60. Even in the I band, 
where the effects of extinction are more moder- 
ate, this diflPerence corresponds to 0.11 mag. As- 
suming that we have chosen the appropriate sub- 
population to match to the Cepheids, the non- 
Gaussian tail creates significant diflFerences be- 
tween the mean and median Ay- For example, 
for the Te > 22000 K sample the mean and me- 
dian Ay are 0.57 and 0.44, respectively. Which 
should one apply? Because of the unknown small 
scale structure of the dust, a photometric correc- 
tion will be statistical and the properties of the 
Ay distribution are critical. As discussed in detail 
by Madore and Freedman (1998), the issue of ex- 
tinction in relation to the distance scale has been 
extensively investigated and the most promising 
techniques are those that sidestep the issue by ei- 
ther using reddening-free quantities or by moving 
as far into the infrared as possible to minimize the 
effect of extinction. Nevertheless, investigators are 
often left with no option but to estimate the ex- 
tinction from "local" stars. We have shown here 
that this approach has possible systematic prob- 
lems. 

Our results suggest another way to mitigate 
the influence of dust. Because the extinction is 
lower for the colder, older stars, using Pop II 
standard candles is analogous to going to the in- 
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frarcd for Pop I standard candles. In a previous 
study (Stanek, Zaritsky, & Harris 1998), we used 
our survey data to investigate the small distance 
modulus, 18.08, obtained by a previous study of 
field red clump stars (Udalski et al. 1998). Using 
the map of HZT we were able to select low ex- 
tinction regions and then correct the photometry 
for the extinction that was present {Aj = 0.31, 
E{B -V) = 0.17). That analysis confirmed the 
small distance modulus determined by Udalski et 
al. . Since that study, the discussion in the lit- 
erature has focussed on the magnitude of the age 
and metallicity corrections to the red clump mag- 
nitude. Corrections derived from stellar models 
(Cole 1998; Girardi et al. 1998) revised the red 
clump distance modulus to 18.3, which is still be- 
low the "standard" value of 18.5 but within the 
uncertainties. The diSiculty with theoretical cor- 
rections is that they are based on models of a com- 
plex phase of stellar evolution. Empirical studies 
of the red clump conclude that the age and metal- 
licities are more moderate (Udalski 1998a, 1998b) 
and that the distance modulus is still ~ 18.1. The 
difficulty with empirical corrections is that they 
depend on the complicated interplay of system- 
atic uncertainties in the various distance estimates 
used to calibrate the age and metallicity depen- 
dence of the red clump luminosity. Our study does 
not address the age and metallicity corrections, 
but rather returns to the issue of the extinction 
correction. 

The HZT extinction map Stanek et al. used is 
based on observations of OB stars and here we 
find that such a map is not an accurate represen- 
tation of the reddening to red clump stars. Our 
low Te extinction map is ideal in correcting for 
the extinction because red clump stars are within 
the Te range included in the construction of this 
map. Our new data suggest that Ay for the field 
red clump stars in regions of low extinction (as se- 
lected by Stanek et al.) is ~ 0.1 and so Aj ~ 0.06, 
which increases the derived distance modulus by 
0.25 mag to 18.32 for the data presented by Stanek 
et al. . If age and metallicity corrections of the 
magnitude suggested by Cole (1998) or Girardi et 
al. (1998) are then applied (~ 0.2 mag), the red 
clump distance modulus can be made as large as 
18.5, in excellent agreement with the "standard" 
value. Whether the latter corrections are neces- 
sary is not addressed by our data and should not 



be inferred from the agreement with the standard 

distance value. 

We do not want to leave the reader with the 
impression that the correction for population- 
dependent extinction removes all of the discrepan- 
cies among LMC distance determinations. There 
remain several argiimcnts for an LMC distance 
modidus that is smaller or larger than the stan- 
dard. Some analyses of RR Lyrae observations 
suggest a small LMC distance modulus (Walker 
1992a; Layden et al. 1996), and an erroneous red- 
dening correction is not responsible for artificially 
decreasing the apparent distance because the red- 
denings for cluster RR Lyrae are derived from 
color-magnitude diagrams. The extinctions to 
these clusters are an accurate measure of the ex- 
tinction toward the RR Lyrae (cf. Walker 1992b 
for an example). Only two of the seven clusters 
used by Walker 1992a have E{B - V) > 0.1, in 
agreement with our generally low extinctions to 
older populations. However not all RR Lyrae 
measurements distances are small, RR Lyrae dis- 
tances using Hipparcos-based calibrations have 
given much larger distance moduli (18.65, Reid 
1997; 18.63 Gratton et al. 1997). Some other 
methods, for example the analysis the SN 1987A 
ring, which again is not affected by the problems 
discussed here, can result in a smaller distance 
modulus (cf. < 18.37 (Gould and Uza 1998)), but 
can also result in larger values (cf. 18.5, Panagia 
et al. 1991; 18.43, Sonneborn et al. 1997). Specif- 
ically for the red clump stars, Udalski (1998b) 
presents data from LMC clusters, for which the 
extinctions are small, and obtains a low value 
of the distance modulus that is in agreement with 
field red clump distance modulus with the original 
extinction correction. The disagreement between 
the field and cluster red clump distances that re- 
sults from our new extinction correction is not 
easily resolved. The clusters have low extinctions, 
so even if the adopted extinction is incorrect, the 
effect on artificially depressing the distance mod- 
ulus is minor. Presumably the cluster stars have 
similar ages and metallicities to the field stars, so 
the resolution will not come from differential age 
and metallicity corrections. 

These unexplained discrepancies still render 
the distance measurement to the LMC uncertain. 
However, we conclude that population-dependent 
extinction is yet another issue that must be ad- 
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dressed as one measures distances and that extinc- 
tion differences between populations can result in 
errors comparable to other, better appreciated, 
sources of systematic error. 

4-6.2. Some Implications for Studies of Stellar 
Populations 

The results presented here also have implica- 
tions for the observation of stellar populations in 
both resolved and unresolved galaxies. The study 
of resolved populations has been advanced in the 
last few years by deep HST color-magnitude dia- 
grams with many thousands of stars, and ongoing 
ground-based surveys that will provide colors and 
magnitudes of millions of stars. Aside from identi- 
fying rare stellar populations that can provide key 
clues to the star formation history of the galaxy, 
the large samples allow precise determinations of 
the apparent positions and dispersions of popu- 
lations, like the red giant branch, which in turn 
can constrain the range of ages and metallicities 
present in the system. The dominant source of 
error in such analyses will no longer be statisti- 
cal but systematic. The extinction correction is a 
principal potential source of such systematic un- 
certainties. We were in fact driven to this analy- 
sis by our inability to fit the main sequence and 
red giant branch simultaneously when applying a 
single distribution of extinctions to a portion of 
our LMC data. Any detailed, precision analysis of 
stellar systems that includes relatively recent star 
formation must account for population-dependent 
extinction. 

The situation is even more complex in unre- 
solved galaxies. As previous studies have stressed 
(cf. Witt & Gordon 1996), the interpretation of 
the measured surface brightnesses is complicated 
by light scattered back into the line-of-sight and 
multiple scatterings. Neither of these issues are 
addressed by the current data on the resolved 
LMC because any scattered light is removed as lo- 
cal background to each star. Ignoring scattering, 
we can estimate the amount of light lost due to 
extinction from both the hot and cool stellar pop- 
ulations in the LMC to determine how an unre- 
solved galaxy with such stellar populations might 
be perceived. In the V band, extinction has re- 
moved 28% of the light from the hot stars but 
only 13% of the light from the cool stars. In the 
B band, extinction removes 33% of the light from 



the hot stars but only 16% of the light from the 
cool stars. Correcting the apparent luminosity of 
the galaxy using the extinction appropriate for OB 
stars (or alternatively from emission line extinc- 
tion measurements (Ha/H/3), which are likely to 
be similar to the OB values), will overestimate the 
galaxy luminosity by ^ 15% and introduce a color 
correction that is about 2 times larger than that 
appropriate for the bulk of the stellar population. 
Although investigators have been aware of such 
issues, the difficulty in dealing with the large un- 
known number of parameters has led to very sim- 
ple treatments. Our discussion here provides an 
indication of the degree of uncertainty that such 
simplifying assumptions may impose. 

5. Conclusions 

We have fit colors derived from stellar models 
to 4-filter photometry of LMC stars to derive ef- 
fective temperatures and visual extinctions. We 
find that there are two populations of stars for 
which we can reliably recover the visual extinc- 
tion, stars with 5500 K < Te < 6500 K and stars 
with Te > 12000 K. Various tests demonstrate 
that we are not strongly sensitive to metallicity 
and that we recover extinctions in agreement with 
previous studies to within a dispersion of ~ 0.1 to 
0.2 mag for individual stars. More precise extinc- 
tion values are derived by combining the results 
from many stars and producing extinction maps 
of the region. Our principal result is that there is 
a significant difference in the extinction values de- 
rived for the two samples of stars, with the hotter 
stars being more highly extincted on average. 

The difference in extinction properties as a 

function of stellar population is not entirely unex- 
pected, but we present direct evidence and quan- 
tify the effect for the LMC. We test for biases by 
comparing to previously published extinction mea- 
surements and to 100/xm images of the region. We 
find that in general the lOO/xm emission features 
trace the extinction morphology, but that there 
are a few exceptions. These few exceptions can be 
traced to the presence or absence of young stars 
that heat the dust. We identify several locations 
where young stars have heated the dust so that 
it emits strongly at 100/xm and locations where 
there are young stars, but apparently no dust be- 
cause there is no 100/xm emission. We conclude 
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that significant 100/ini flux correlates with young 
stars, but that young stars do not necessarily im- 
ply significant 100/xm flux and that 100/xm flux 
should only be used as a star formation tracer with 
caution. This is in agreement with the results of 
Walterbos and Schwering (1987) for M 31. 

We presented a simple model for the distribu- 
tion of stars and dust along the line of sight. Pre- 
vious results (HZT) demonstrate that the scale- 
height of the dust is about twice that of the OB 
stars. Because the low Te stars have different ex- 
tinction properties, we hypothesize that the scale- 
height of this populations is ^ than that of the 
dust. Such a model implies that half of the low Te 
stars are completely in front of the LMC dusty disk 
and that half are extincted by the extinction at 
that location in the disk (given by the high Te ex- 
tinction map). We then compare the expectation 
from this model with the spatial correlation be- 
tween low and high Te extinction measurements. 
We find that this model fits the data remarkably 
well and provides a simple explanation for the dif- 
ferent extinction properties of the two samples. 

We explore some ramifications of the population- 
dependent extinction. In particular, we discuss its 
effect on distance measurements of the LMC. The 
expected behavior is in the correct sense to help 
resolve the discrepancy between distance mod- 
uli found using various stellar populations (cf. 
RR Lyrae, red clump stars, and Cepheids). We 
present how the distance measured using field red 
clump stars may have been biased downward by 
using reddening measurements obtained from OB 
stars, rather than from red clump stars, in the 
same region. Taking the mean extinction obtained 
from the cooler stars, which are directly compa- 
rable to the red clump stars in temperature, the 
distance modulus inferred is 0.2 mag larger than 
previously derived, in better agreement with val- 
ues obtained from observations of Cepheids. How- 
ever, the Ccpheid distance measurements may 
also be affected by population-dependent extinc- 
tion variations because young stars are subject 
to greater and more variable extinction, making 
determination of the Ay distribution, rather than 
a mean or median, key in properly accounting for 
the effects of extinction. We discuss how observa- 
tions of galaxies in general may be affected by the 
population-dependent reddening. Treating this 
section of the LMC as a complete independent 



galaxy, we derive that nearly twice as much of the 
light from the hotter stars is removed by extinc- 
tion than of the cooler stars. This effect results in 
an effective extinction curve steeper than the true 
curve. 

We conclude by stressing that dust affects dif- 
ferent stellar populations differently. Simple mod- 
els of obscuring planes or well-mixed stars and 
dust do not accurately represent the physical sit- 
uation. Whether this gravely affects our interpre- 
tation of observations is still not fully understood. 
We illustrate that in certain cases ignoring the 
results presented here could have serious conse- 
quences. We continue to reduce the data obtained 
for the LMC and SMC as part of our survey and 
soon should be able to examine the distribution 
of dust on even larger scales. Then we can ad- 
dress whether the lack of dust in the center of the 
starforming region in LMC-4 is common for bub- 
bles and superbubbles, how the dust is distributed 
though the disk of the LMC, how it is distributed 
in the bar, and whether similar distributions are 
observed in the SMC. These data provide a unique 
opportunity to explore the relative distributions of 
stars and dust within galaxies and are a necessary 
first step in properly interpreting the highly popu- 
lated color magnitude diagrams that are becoming 
available for the Clouds. 
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Fig. 1. — The stellar density plot (for stars with V < 21) for the region in the LMC with reduced photometry 
from the Magellanic Cloud Photometric Survey. The central coordinates are roughly a = 5^20^ and 
5 = -66°48', the image is - 4° wide by 2.7° tall, with North at the top and East to the left. Each "pixel" 
in this plot corresponds to 21". 
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Fig. 2. — Comparison of photometrically derived effective temperature with spectroscopically derived effec- 
tive temperature. Each panel depicts the results for standard stars from the application of an algorithm that 
matches stellar models to data. The results in the upper left panel were obtained with neither a correction 
for extinction or a photometric zero point shift between the standard stars and the models. The results 
in the upper right panel were obtained from models that include extinction as a fitting parameter, but do 
not include a zero point photometric shift. The results in the lower left panel are from models that include 
extinction and photometric shifts. The results in the lower right are from the application of models that 
include extinction and the photometric shifts to an independent set of standard stars. In the upper left of 
each panel we provide the fractional shift between the two temperatures and scatter from the 1:1 line. 
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Fig. 3. — Recovering a supplemental, artificial line-of-sight extinction. The top three panels contain the 
recovered extinction {Ay) as a function of recovered effective temperature (a running median in Te for 11 
stars) drawn as a solid line. To the existing data, wc added a visual extinction of 0.2 mag to generate the 
upper panel, 0.4 for the second panel, 0.8 for the third panel. Each panel has a dotted line that corresponds 
to the level of supplemental, artificial extinction. The lowest panel shows the extinction recovered from the 
simulation with a supplemental extinction of 0.8 (same as the third panel) but the data are plotted versus 
spectroscopic Te- 
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Fig. 4. — The distribution of visual extinctions for the original sample of low Te stars (solid line) and for 
the corresponding sample with larger photometric errors (0.2 < au < 0.4; dotted line). The hashed areas 

represent the stars for which the derived Ay = 0. The two histograms have been arbitrarily normalized to 
match at Ay = 1. The original sample is 3.4 times larger than the high au sample. 
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Fig. 5.— Extinction maps from stars with 12000 K < Te < 45000 K (top) and stars with 5500 K < Te < 
6500 K (bottom) for the region showed in Figure 1. The lower image has been smoothed with a Gaussian of 
1.5 map pixel dispersion to clarify any structures. The brightness of the images are identically scaled (from 
Av = (white) to Ay = 0.75 (black)). The small vertical bar to the right of the images marks the position 
of a scan known to have non-photometric /-band data. 
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Fig. 6. — The dependence between metallicity and the derived extinction and effective temperature. The 
mean and dispersion of the difference between recovered quantities for various models are plotted. The 
comparison of results using standard models ([m/H] = —0.5) to solar metallicity models is shown in the left 
panels. The comparison to tenth-solar metallicity models in the right panels. The open circles at Te ^ 6000 
K represent the median values in those bins (see text for discussion). Errorbars are the dispersion of the 
mean. 
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Fig. 7. A comparison of extinction values derived for stars in common with the HZT study. We compare 
our values to the values they derived using the classic UBV reddening-free method. The line indicates the 
one-to-one correspondence and is not a fit to the data. 
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Fig. 8. — A comparison of the reddening measured by OSG to Galactic stars and the foreground reddenings 
implied by our map at the positions of the stars. Although the mean and dispersion of the two populations 
(i.e. the global statistical quantities) agree well, there is no discernible coherence in the point-by-point 
comparison. 
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Fig. 9. A comparison of the distribution of extinction values derived for the subset of our cool stellar 
population, 5500 K< Te < 6000 K (solid line) that is interpreted to be predominantly Galactic foreground 
stars and for Galactic foreground stars (dotted line) from OGS (1995) over their entire LMC region. The 
hashed box represents the number of stars with Ay = in OTir analysis {Ay < is not allowed in our 
models). The histograms are arbitrarily normalized. In the inset we present the Ay distribution for both 
our Galactic foreground stars (thick line) and LMC population (thin line). 
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Fig. 10. A comparison of the foreground reddenings determined from a map construeted from the OSG 
data and our extinction map constructed from all of the cool stars (LMC and Galactic). Pixel values are 
binned according to the OGS extinction value and averages of the pixels from our map are plotted (errorbars 
are the dispersions of the mean). The correlation is evident, indicating that structures seen in the Galactic 
foreground map are also evident in the LMC data, despite that lack of an apparent correlation in the 
extinctions of individual stars (Figure 7). 
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Fig. 11. — A comparison of the Ay map derived from high Te stars (smoothed to match the IRAS spatial 
resolution), an IRAS lOO/xm image (log scaled), and a stellar density image of main sequence stars with 
V < 16.5 and B — V < 0.5. The vertical bar at the right of the upper panel indicates the scan that has 
non-photometric /-band data. 
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Fig. 12. — A comparison of the distribution of Ay's derived for the hot (Tg > 12000 K; thick sohd hne) and 
cool (5500 <Te < 6500 K; thin soUd line) stars. The dashed areas at Ay < indicate the number of stars 

with Av = 0, where wc have set the height of the bin equal to the number of stars with Ay = 0.01 and the 
area to correctly represent the number of stars. The two histograms have been arbitrarily normalized. 
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Fig. 13. — Comparison of extinction values for subsets of the high Te population. The lower two panels 
compare the values obtained by HZT using the reddening-free method. 
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Fig. 14. A comparison of the extinction derived for LMC stars with 12000 K < Te < 45000 K and stars 
with 5500 K < Te < 6500 K. The data are map pixel values binned so that there are at least 200 pixels per 
bin. Errorbars represent standard deviation of the mean. The lines represent models as described in text. 
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